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Abstract 


By examining the exact analytic solution of a kinetic model of colllslon- 
al Interaction of Ionospheric Ions with atmospheric neutrals In the Bhatnagar-* 
Gross~Krook approximation, we show that the onset of Intense auroral electric 
fields in the tppslde Ionosphere can produce the following kinetic effects: 
(1) heat the bulk Ionospheric lone to ~ 2 eV, thus driving them up to higher 
altitudes where they can be subjected to colllslonless plasma processes; (2) 
produce a non-Maxwelllan superthermal tall In the distribution function; and 
(3) cause the Ion distribution function to be anisotropic with respect to the 
magnetic field with the perpendicular average thermal energy exceeding the 
parallel thermal energy. 
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1. INTRODUCTION 


Ground-based and satellite observations in the latter part of the past 
decade have brought about a major advance in the understanding of auroral arc 
formation processes: the electrodynamic interaction between the hot magneto- 
spheric plasma and the cold ionospheric plasma seems to play a central role 
te.g., reviews by Akasofu, 1981; Mozer, 1981; Kan and Lee» 1981; Chiu et al.| 
19811. Although a major consequence of this electrodynamic interaction is the 
production of a component of the electric field parallel to the magnetic field 
for the acceleration of electrons in discrete arc formation, a second, and 
possibly more far-reaching, consequence of the Interaction is the finding that 
the ionosphere is a significant source of plasmas in the magnetosphere as 
ionospheric ions are accelerated upwards by auroral electric fields related to 
substorms te.g*, Shelley et al., 1976; Mizera and Fennell, 1977; Richardson et 
al., 1981]. This is in addition to the polar wind [Banks and Holzer, 1968; 
Banks, 1979] which Is a significant steady mechanism for transporting iono- 
spheric plasma into the magnetosphere. Actually, the idea that the ionosphere 
is an active participant in the magnetospherlc response to solar-terrestrial 
activity has had a fairly long history [e.g., Dungey, 1961; Axford and Hines, 
1961; Vasyliunas, 1970]. These authors pointed out the importance of the 
ionosphere in providing the appropriate Pedersen conductlvi/Ay to close the 
magnetospherlc convection circuit - a role which is crucial in modern theories 
of magnetosphere-ionosphere coupling [e.g., Chiu and Cornwall, 1980; Kan and 
Lee, I960]. In connection with this role of the ionosphere in limiting mag- 
netospherlc convection. Joule heating of the ionosphere [Walbrldge, 1967; 
Redder and Banks, 1972] and the neutral atmosphere [e.g., Chlng and Chiu, 
1973; Straus and Schulz, 1976] have been considered not only as a thermal 
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onergy source but also as a source to drive ionospheric motions which have an 
indirect influence upon the raagnetospheric convection flow# These studies do 
not address the question of the ionosphere as a substorm’'related source of 
magnetoBpheric plasma* In this paper we shall attempt to consider the kinetic 
properties of the ionospheric plasma as it responds to the onset of an 
enhanced auroral electric field, and as a source of raagnetospheric plasma. 

The observational picture of the kinetic plopertles of such magneto" 
spheric plasmas of ionospheric origin is far from complete; thus, any serious 
theoretical effort at present must be in the category of "base building." 
Upward acceleration of ions in auroral electric and magnetic fields is probab- 
ly not difficult to understand since both the parallel electric field in 
Inverted-V structures and the divergence of the magnetic flux tube favor 
adiabatic upward acceleration of ionospheric ions such as 0*^. Nonadinbatic 
features of auroral plasmas of ionospheric origin (such as heating and genera- 
tion of Buperthermal populations of Ion beams and conics) are an entirely 
different matter* For Ionospheric ions, the observations are particularly 
intriguing since not only are these ions somehow energized to superthermal 
energies in directions parallel [e.g., Richardson et al., 1981] and perpen- 
dicular [e.g., Klumpar, 1979] to the magnetic field, but the processes seem to 
operate over very wide ranges in energy (-* 6 eV - iO keV) and In altitude (500 
km - 8000 km), and over a wide distribution of local times [Gorney et al., 
1981 j. Such preliminary observational results clearly indicate the direction 
of present and future theoretical studies of auroral Ionospheric ions: how 
are superthermal and nnisotroplc ion populations formed /and hovj are ions 
energized over some four to five orders of magnitude in energy? 

Quite posslb3,y, the answer to these questions may be in-sltu wave- 
particle interactions le.g*, Un^wtrup et al., 1979; Okuda and Ashour-Abdalla , 
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1981]; perhaps future quantitative simulationB with more realistic conditions 
will answer the question of formation of beams (parallel energization) and 
conics (perpendicular energization) in the observed energler covering the 
ranges of 6 eV to > 10 keV. Whether such theories are realistic or not is 
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outside the concern and scope of this paper. Rather, we are more interested 
in the origin of such superthermal populations at the lowest energy ranges (< 
6 eV). The question of how special nonadlabatlc features are formed out of 
the cold ionospheric ion population is quite puzzling when one examines the 
conditions of the ionosphere. Flvrst, hlgh-latltude ionospheric temperatures 
at F-reglon heights, as measured by the S3-3 satellite, are generally < 2500“ 
K ('^0.2 eV) [Rich et al., 1979], far less in energy than the superthermal 
fluxes observed. Second, Kindel and Kennel [1971] concluded that the unstable 
regions of electrostatic ion cyclotron waves were above the F*~maximum and 
generally in the far topside ionosphere above IQOQ km; this makes the relative 
abundance of 0*^ in such nonadlabatlc populations (conics and beams) 
[Ghlelmetti et al., 1978] even more puzzling since must be driven up to 
such altitudes by some preheating process so that 0 can at times be the 
dominant ion at altitudes > 1000 km. Note that in this situation, the abun- 
dance of O'*’ at altitudes > 1000 km cannot be attributed to escape of H**" to 
higher altitudes, as 0"*’ is also observed to be dominant at the equatorial 
regions during storm time [Balslgdr et al., 1980]. Could it be possible for 
some kind of pre-heating process to operate in the ionosphere to drive up an 
intense superthermal (~ several eV) population of ionospheric ions into 
regions where ion cyclotron waves can act to energize these ions up to tens of 
keV? Recent observations [Lockwood and Tltherldge, 1981] support such 
hypotheses. 
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As we have discussed hithecto, the probable occurrence of ion Joule 
heating in the ionosphere by convection electric fields mapped down to iono- 
spheric heights tf odder and Banks, 1972) Ims been well accepted. Since such 
calculations use a fluid approach, which automatically assumes a Maxwellian 
form for the distribution function, they 4^nnot tell us about kinetic features 
of the ion distribution such as pitch angle anisotropy and superthermal popu- 
lations without going into extremely complex calculations with higher moments 
[Schunk, 1975). In view of the necessity to understand how magnetospheric 
ions may originate from the ionosphere and in view of the predominance of the 
steady convection electric field as an ionospheric heat source, we are 
prompted to ask if the auroral electric field may not be the source of pre- 
heating which provides the topside with superthermal ions. Since the observed 
nonadlabatic features are ki netic in character, we are then driven to co ns ider 
the kinetic response of ionospheric ions to the onset of steady auroral elec- 
tric fields in a ‘simple colllsional kinetic plasma model - the Krook model. 
Analytically soluble models always oversimplify; thus our main purpose is to 
gain some Insight Into kinetic properties of the ion distribution under 
auroral ionospheric conditions rather than to attempt to “explain" the 
observed nonadlabatic properties alluded to above. Our work, simple though it 
may be, differs from previous work on ionospheric heating using a fluid 
approach [Fedder and Banks, 1972) in that we are required to examine the non- 
Maxwelllan features of the distribution function rather than to assume 
Maxwellian distributions at all times, as is done in fluid models. 
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II. formulation and preliminary considerations 

Consider a unifoiro auroral eleetrlc field imposed upon a horizontally 
uniform ionosphere consisting of ions, electrons and neutrals whose distribu- 
tion functions are respectively J fj^, fg and f^^. In the auroral ionospheric 
region of interest, where we assume the magnetic field to be vertical and 
uniform - Bz), these thr/ee species interact through collisions and the 
total self -consistent electric field The Interactions of these species are 
symbolically formulated in terras of three coupled Boltzmann aquations for the 
distribution function fj^, k-(i,e,n){ 


L 


£ 

k \ 


E 

i^k 




(i) 


where 




io the Boltzmann operator for species k of charge qj^e and mass m;^, 


= . v+ 5 D/Dt 


( 2 ) 


The force consists of a combination of electric, magnetic and gravitational 
components. 


tk"\8’*'®*lk^*^® qj^(v X B z)/c 


( 3 ) 


is the appropriate binary collision operator. The self-consistent elec- 
tric field must satisfy Poisson's equatloa, 


7 • t ■ 4ne / d^v (f. - f ) 

X 6 


( 4 ) 
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It ii weU-knovm Chat (1) - <4) are very difficult to eolva even with the 
simpleat assumptions for K# In the next section we shall show that a trun*" 
cated form of these equations with a Krook-type collision term can be solved 
to give some insight into the kinetic properties of ionospheric ions under the 
Influence of a strong Ej^* 

The highly symbolic formulation (1) - (4) nonetheless allows us to 
discuss certain asymptotic constraints relevant to the problem. If we ignore 
the fluctuation part of IS <the auroral DC is usually observed to be much 
larger than the AC electric field), the drivers of non-adlabatic effects are 

A 

the non-linear collision operators which allow energy to be exchanged 
among the three species. Thus, given sufficient time and no loss of particle 
species from the system, it is believed that the colllslonal evolution of the 
system tends asymptotically toward an isotropic Maxwellian form for the 
distributions of all the species, although the equilibrium thermal energy of 
each species need not be the same. This means that the fluid moments 
approach, with the underlying assumption of isotropic Maxwellian distribu- 
tions, is valid at times long after onset of for ions the time scale is 
likely to be set by the ion-neutral collision time. In the E-region, where 
ionic Joule heating effects have been thoroughly studied [Fedder and hanks, 
1972] using the fluid equations, the equivalent time scale as determined by 
the ion-neutral collision time is much less than one second but the neutral 
response time is ~ 1/2 hour; therefore, we do not expect kinetic features such 
as superthermy and pitch-angle anisotropy to persist for much longer than the 
ion-neutral collision time after onset. The situation for the F-reglon and 
the topside (> 350 km) is quite different. In this regime, the ion-neutral 
collision time (> 10 sec) is much longer than the ion gyration period (~ 1/30 
sec); therefore, the expected evolution, to Maxwellian isotropy takes place 


6 


over many gyration periods, perraicttog kinetic features to persist for tens of 

seconds or minutes. Can a strong drive superthermy and pitch-angle 
* 

anisotropy in the ionospheric ion population during the evolution period 
between Ej^ onset and the asymptotic Maxwellian state? It is the intent of 
this paper to demonstrate that this question con be answered in the affirm- 
ative for the Krook model of weakly ionized plasma kinetics. 

Before we proceed with the solution of the tonic segment of a Krook-model 
plasma in uniform electric and magnetic fields, it may be convenient to give 
an elementary discussion of the origin of the expected kinetic features of 
superthermy and pitch-angle anisotropy. An imposed Ej, causes ions to drift 
and collide with neutral atoms of the thermosphere which act to deflect the 
uniform drift motion into somewhat random motion. In the Initial stages (far 
from the asymptotic state) the rate of change of ion kinetic energy is not 
isotropic because the mobilities (iC|, of ions in the F-region are 
different for motion parallel and perpendicular to the magnetic field. Tl\e 
perpendicular energy changes according to [Alfven and Falthammar, 1963] 


. ‘"'ll/'*" ■ ''i - “l '’in 


(5) 


where [Rlshbeth and Gatriott, 1969] and « 2/3 for the 
F-reglon. In (5) and Wj^ refer to mass and energy of the neutrals; is 
the ion-neutral collision frequency and ft is the ion gyrofrequency. The first 
term on the right hand side of (5) is due to the driving of the electric field 
and the second term is due to energy loss to the neutrals. The same amount as 
the energy loss in (5) appears in the corresponding equation for as neutral 
energy gain. The simultaneous solution of (5) together with the equivalent 
equation for does not concern us here. The point that the ion population 
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may be anisotropic soon after Ej^ onset is made by notlitg that the parallel ion 
energy is unlikely to Increase until sufficient perpendicular drift energy has 
been scattered through large angles into the parallel direction - probably 
after may coUialoa times; whereas, by virtue of (5), Inoreaaea In one 
or a few collision times ♦ 

Obviously, the ideas d^scussed in this section are qualitative. For the 
rest of this paper, we shall consider the results of a simple solution of an 
ion component of the Krook model for weakly ionized plasma in steady electric 
and magnetic fields. 


III. MODEL OF KINETIC RESPONSE 


The TuU treatment of the kinetic response of a gravitationally- 
stratified weakl;,”*ioni 2 ed and inhomogeneously-magnetized plasma (the F-region 
ionosphere) to the onset of a strong electric field, as formally presented in 
the previous Section, is exceedingly difficult to solve. It is the purpose of 
this Section to initiate a kinetic analysis of this transition region between 
the collision-dominated E-reglon and the co'' Hsionless magnetosphere with on 
oversimplified but exactly soluble version of (1) - (4). The simplifications 
Introduced are; 

1. The neutral component is assumed to remain in static equilibrium. 
Thus, the neutrals act as scatterers of ion motion but do not pick up any 
energy in the process. This is a reasonable approximation if the ionic con- 
centration is sufficiently low or if the average neutral mass is high. 
Neither condition is strictly valid for the F-reglon. Redder and Banks [1972] 
showed that in a fluid model the motion of the neutrals is an Important 
determinant of the ion heating on time scale greater than about one half hour 
after electric field onset; since we are dealing with ion response at early 
times (several collision times after onset), the assufiiption of Immobile 
neutrals is approximately valid. ‘ We hope to relax this restriction in our 
nes:t stage of kinetic model development. 

2. We Ignore the colllslonal Influence of electrons upon the dynamic 
response of ions to steady electric field onset. This is not a bad 
approximation for early times because ion-electron collisions do not change 
the ion energy by very much. 
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3. Wq Ignore the fluctuating (AC) electric field in (^)# The AC field 
driven by auroral electrona may be very important in ion heating, especially 
in the form of resonance with electrostatic ion cyclotron waves. This effect 
has been pointed out by others (e.g», Okuda and Ashour-Abdalla, 1981). Here 
we are not interested in cyclotron heating but, as we have stated, concentrate 
primarily on how the bulk of auroral ionospheric ions at ~ 0.2 eV may be 
heated to several eV on the topside. A complete thoory of ion conics and 
beams cannot, however, ignore the effects of the AG electric field. 

4. The collision operator for the ion component of (1) is approximated 
by the simple Krook model [Bhatnagar et al., 1954], which was proposed speci- 
fically to study the approach to equilibrium of weakly-ionized collisional 
plasmas. 

*^in " ^ ^ (x,v,t) + n (x,t) f^ (x,v,t)] (b) 

where V is the ion-neutral collision frequency which thermalizes the ion 
distribution function f towards the aseuaied isotropic Maxwellian form fg 
[Bhatnagar et al., 1954], 

fg s [ra/2r T(x,t)]^'^^ exp - [m/2T(x,t)] (v - u(x,t)]^ (7) 

Note that / d'^v fg - 1. In (6) and (7) n, u and T are the density, flow and 
energy moments: 


n 5 / d^v f 

(8) 

Us / d\ V f/n 

(9) 
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T « / d^v m[v - u]^ f/3n (10) 

In order for the Krook model to conserve the above three quantities, (6) - 
(10) have to be self-consistently included in the solution of (1) [Bhatnagar 
et al., 1954]. The ion-neutral collision frequency v can depend on n, but we 
shall assume v to be a given function of sa only. By using the Krook model, it 
is assumed that v does not depend on the velocity v because the Krook moi, 1-^.1 
does not conserve number and energy if v v (v, x). 

5. The plasma is assumed uniform in the horizontal direction and g * 0 
in (3). Obviously, this simplification Ignores auroral spatial scales, but 
little progress can be made otherwise. The distribution function f can thus 
depend only on z, the vertical coordinate. 

6. The driving DC field ^ is assumed uniform for t > 0 

I - Ij^(t) - 0(t) X (U) 

Since V • ]£ » 0,’ (4) implies / d\ f " ®0 ionospheric electron 
density - assumed constant. Thus, seif-consistency of the electric field (4) 
and assumption (11) implies n(x,t)’" uq. This constraint can be verified from 
our solutions. 

Our model of kinetic response is defined by (1) - (4) and (6) - (11). 
With the minor exception of restrictions on t and v, this model is essentially 
the extension of the Krook model to the case of collisional plasma in uniform 
electric and magnetic fields. We seek solutions of this model as an initial 
value problem in which ^j^(t) causes f to evolve with time from an initial 
isotropic Maxwellian state 
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f(x, V, t 0) « exp - (m/2TQj 


( 12 ) 


where Tq Is the cold Ionospheric temperature (~ 2500“ K)j T (x,0) »« Tq, 
u(x» 0) « 0. 

Since the assumed force f is divergenceless in space, we seek a solution 
of our model by a phase-space transformation (x, v, t) (x’, v', t) so that 
D/Dt in (2) is transformed into 3/9t. 'ihis is accomplished by the following 
transformation between the components, labeled (1, 2, 3), of the above 
vectors: 


^3 “ ^3 


[''11 


~cos nt - sin nt 


V, 


- sin nt 


m 



X 

•f (cEj^/B) 


/2j 


sin nt cos nt 

*■ 


72. 


cos nt - 1 


e(t) 


(13) 

(14) 

(15) 


*1 

m 

I 

X 

1 

- Iq 

v}(t) 





. - 


e(T) 


(16) 


Note that the coordinates (x', v! ) are the time-reversed evolution of the 
coordinates of an ion under and Bz starting at (x, v) at t ■ 0; thus 


v‘(x. o)].rn 
X’(x, V, 0)J [x] 


(17) 


Further,' since the transformation (14) between velocity spaces entails a 
simple rotation, the volume element is invariant: d^v ■ d^v' . Tl»e transla- 

tion in (14) amounts to a shift of origin in velocity spacevi Under this 
transformation, the Krook equation Ifr/iim (1), (2) and (6)] 
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( 18 ) 


becoroQB (see Appendix A) 


Df _ y 
Ut *^ln 


3f(^,y',t)/3t - v[z(^,^,t)]{nQ.?Q[x(x’,v',t), v(^,v',t),t] - £(x' ,v’ ,t)} . (19) 


In (19), for the sake of explicitness, we have written x, v as functions of 




(x’, v’, t). Since we shall be working with the (x‘, v*, t) coordinates we 


shall hereafter use x and v to denote functions of (x', v’, t) specified by 
the inverse transformation of (13) " (16). 

The solution to (19), with the initial conditions (12), is 




f(x',v’,t) - f(x\v’,0) e”''*^ + £q[x(x’ ,v' ,t), v(x',v't),t] (20) 

U 0 




where v E v[z(x’ ,v' ,t)] in general. Note that this integral representation 
solves (19) generally once f(x’,v',0) and fQ(x,v,r) are given. General inte- 
gral representations for u(x',t) and T(x',t) are obtained by substitution of 

I 

(20) into (9) and (10), Because of assumption 6 (i.e., * 0 and (4)), the 

possible forms chosen for fQ are constrained by 


n(x',t) - / d^v’ f(x',v',t) ■ Uq 


We show in Appendix B that our choices (7) and (12) satisfy this constraint 
(at least for the simple model given below). The integral representation (20) 
shows explicitly that, as the consequence of collision with neutrals, the 
initial distribution f(x' ,v’ ,0) disappears in a collision time (1/v) while a 
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new convoluted distribution takes Its place* 


We have. In principle, solved the problem with (20) and (21) j however, 
devising a tractable model so that (20) can be evaluated explicitly is another 
matter, since the complexity of the problem is now hidden in the transformed 
fQ and V. We have examined a number of models and the most tractable requires 
that V « constant, in evaluating the moments u and T, (9) and (10), which are 
crucial in the solution (20) because fQ depends on u and T. For the rest of 
this paper, we shall restrict our discussion to the "simple" model: v » 
constant* 

In this "simple" model, the explicit solution is reduced to solutions for 
u and T of the following set of coupled integral equations derived by substi" 
tutlon of (20) into (9) and (10): 


u^ 


u(x’t) e ' “ / d v’ V f (x' ,v’ ,0)/nQ 


+ V / dt e'^^ / d^v' V ffs(x,v,T) 


( 22 ) 


‘ [3T(x' ,t)/m] ■ / d\’ [v f(x*,v’,0) 


+ \) j *dT e'*^ / d\' (v - u]^ f^ (x,v,t) 
0 ^ 


(23) 


As we have stated earlier, the dependence of v and x on (x’,v’,t) are sup- 
pressed in (22) and (23) for brevity. To reduce these integral equations into 
tractable form we observe that the velocity space integrals can be performed 
by applying the relations (13) - (16) to (7) and use the identity 


*.2 


(v - u) - [vj - Ujj sin fit “ Uj^ cos ftt + U 2 sin ftt] 
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+ [v^ + Ujj (1 - C08 nt) 


Uj sin nt 


Ug cos nt]^ + [vj - \x^]^ (2A) 


whero uq ■ cEj^/B Is the drift speed. Carrying out the velocity space inte- 
grals of (22) and (23) is somewhat tedious; for purposes of lllustra'iCion, we 
will exhibit only the procedure for ( 22 ) j U 3 •• 0 and 


“1 

t \I*I* 

cos (lx 

1 

•“ (eE,/m) / dT e 



0 

-sin At 


. ( fsln (It 

' I [cos Ot 


cos fit"] rnl -vt 
sin ntj |_vj ® 



(25) 


It is Interesting to note that (22) is decoupled from (23) because it turns 
out that / d v’ V (x,v,t) » u (x*,t), giving the simpler results (25). At 
t««* , Uj and U 2 approach the well-known Pedersen and Hall drifts respectively, 
as we would expect. Using (24) and (25), the Integral equation (23) for the 
thermal en^^rgy 3X (x',t)/m can be written explicitly as 


(3T/ra) - (3 T^/m) + 2 {(1 - e"^''*')/2 - (1 - e”^''*')/3}/(v^ + fl^') 

+ 2up V n {-e"'’^ slnftt 4- cosftt) + (2v^+ n^) slnAt]/ 

(4v^ + n^)}/(v^ n^) (26) 

The thermal energy of the ions is the sum of the initial thermal energy 3Tq/hi, 
a positive monotonic contribution due to colllslonal conversion of the elec- 
tric drift and a term oscillating at the ion cyclotron frequency. The oscil- 
latory term is unimportant in the F-reglon because it averages to order 
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(v/fl)^ which l8 much less than I* As t 

(3T/m) (3 Tq/hi) + y n^/(v^ + (27) 

2 

Thus, about two thirds of the drift kinetic energy mUjj/2 goes into therrjal 
energy. The exact amount of energy conversion depends on model assumptions; 
we regard the result of this model only as a guide to ion heating in the F- 
region. According to (27) then, the thermal energy increase of the bulk of 
topside F-region ions under the influence of the auroral electrostatic field 
is roughly 

(T “ Tq) ^ mc^ (Ej^/B)^ ^ 8 eJ X lO"^ eV (28) 

where is in units of mV/m. Thus, if Ej^ «>- 500 mV/m [Mozer, 1981], F-region 
oxygen thermal energy will be ~ 2 eV - sufficient to energize bulk oxygen 
expansion into the topside auroral ionosphere [Lockwood and Titherldge, 1981], 
but not sufficient to cause O"*" escape which requires > 10 eV. Note that this 
meclianlsm is important only on auroral field lines where Ej^ is large; if we 
consider the large-scale convection field (E^ < 100 mV/m) the thermal energy 
increase is only < 0.1 eV (< 10^ K). This result is in agreement with the 

fluid calculations of Fedder and Banks [1972]. Further, this mechanism favors 
thermal energization of oxygen over hydrogen as observed by Ghielmetti et al. 
[1978]. 


IV. SUPERTHERMY AND ANISOTROPY 

We presented the formulation and solution of the kinetic model in the 
previous section with a discussion of thermal energy increase of F-region 
ions. Since T is the thermal erergy moment of the entire distribution func- 
tion, as defined in (10), the above discussion is not applicable to considera- 
tion of non-Maxwellian features of the distribution function. These are 
discussed in this section. 

To show that the solution (20) contains a superthermal ‘'tail" which is 

non-Maxwellian, we need to express f(v’,t) in an explicit form. This is very 

difficult, since the time convolution integral of the transformed fp is very 

complicated because u and T are complicated functions of time. Instead of 

giving a detailed numerical study, we shall show that a super thermal non- 

Maxwellian "tail" appears at early times: vt « 1 and fit » 1. Since Ot » 

1, the oscillatory terms in u and T can be Ignored because they average to 
2 

order (v/fl) or smaller. Thus, to order (v/fl), (25) and (26) yield 

tij « c(Ej^/m)(v/n) 5^2““ c(Ej/m)j u^ « 0 (29) 

T/m « Tq/bi + ^ Up V^t^; Up «« cE^^/B (30) 

12 

Defining T^/m = -j Up, we can write fQ (v*,t) as 

fp « (m/2TTTp)^/^. (1 + Tg V^t^/Tp)"^/^ • “ I (v»-u)^/(l + Tg v^t^/Tp) (31) 

For (T^/Tq) (vt) « 1, which la somewhat more restrictive than vt « 1, we can 

3 

expand £q above and write an approximate expression of f to order (vt) , 
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( 32 ) 


f « Kt-OKe"''*' + /^dT ^ ^VTj)(VT)^fm(v'4)^- 3 Tq1) 

- f(t«0)U (vt)^ (Tg/T^)(m(^-u)^- 3XgJ + 0[(vt)^J) 

where f(t-O) 1 b the Maxwellian (12). The factor in curly brackets in (32) 
Indicates clearly that £ is a non-Maxwelllan with a superthermal tall in the 
supertherpial regions of velocity space where m(v'“u) > 3 Tq. Note that par“ 

tides in the m(v'-u)^ < 3Tg regions have been shifted to the superthermal 
region. 

A second important feature of the early“tlrae approximation (32) to the 
distribution function is that the superthermal factor (curly bracket) Is 
anisotropic because of (29). The anisotropic factor (v'-u) is thus roughly 

- u)* « tvj - c(Sj^/B)(w/n)J^ + [vj + o(Ej^/b)1“ + (33) 

Because v/n « 1 in the F-reglon, we expect the superthermal part of the 
population in velocity space to be enhanced in the x 1^ drift direction; 

i.e., the 2-direction. More discussion on the anisotropy of the superthermal 

I 

population will be given in the next section. 

The property of anisotropy is not limited to the superthermal part of the 
population. Indeed » we now show that the thermal energy of the bulk of the 
distribution is also anisotropic with respect to the magnetic field. We split 
the thermal energy momenti (10) or (23) » Into perpendicular and parallel 
components : 


"i = ! - “i>^ ‘/"o 


(34) 
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W| i / d^v* (v^ - f/nQ 


(35) 


W i 3T/m " + Wj 


(3d) 


where W is given explicitly by (26). A tedious calculation gives 

2 

W| - (2 “ e'^‘^KTQ/m) ([n'Vs + 2vV/(4v^4(^^)] « 

V 4(2 

4 2v^ (l-cosnt)] e~'*’^ + [n^ - (n cosnc n* 2v sinfit) 2v^(2/(4v^4f2^)] e 


- (n^/3) e“^'’^> 


(37) 


from W and *' 


t« »i c /»s<i\ irtu_ asA"* 4...* fnc\ I /o'TN j..— 

W| we can oocain irota v-»w/» *ims c;f^%pj.(s?i.4.u/ ux v*w/ auu \.~ir / uuc 

to the rapidly oscillating terms is non*-esBenclal. We can set the sinflt and 

(1 - cos(2t) terms to zero by averaging over a time long compared to cyclotron 

time but ^hort compared to collision time and denote the time~averaged thermal 

moments with an overbar. Further* with the approximations v « n and 3TQ/m 


« Up I we obtain- the thermal anisotropy ratio 


Wjl/( 2 W|) » 1 4.|. e’'^*^/(l “ e”'**^) J £9^0 (38) 

where the factor of 2 associated with W| compensates for the two degrees of 
freedom for the perpendicular component. Note that (38) is valid for t > 1/(2 
from onset because we have taken cyclotron averages. For t-0, the exact 
expressions yield Wj^/(2Wj) " 1, as we would expect. From (38)* we note that 
as t +■ « the anisotropy again disappears. During the transition phase 
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(vt ~ 1), th© thermal anisotropy can be quite large 3 or 4). In the 
topside F-regioni the collision times can be tens o£ seconds to minutes » which 
are of the order of the expected lifetime of auroral acceleration potential 
structures [Chiu and Schulc» 1978; Chiu and Cornwalli 1980]* Thus^ we would 
expect thermal ions in the topside auroral F-region to show preferential bulk 
heating in the perpendicular direction* 
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V, DISCUSSION AND CONCLUSIONS 

BecfflUio of the necessity to make simplifying essuroptions in order to make 
the model tractable, we do not claim any direct relevance of our model results 
to observations except in a qualitative sense t However, oven at a qualitative 
level, it is perhaps worthwhile to discuss the implications of our results, 
which are listed as follows: 

1. In the weakly collislonal regime of the topside F»-reglon (v « fl), a 
local steady (100 - 1000 mV/m) can heat the ions to a temperature of 
several electron volts, thus driving the bulk of thermal O'^ to high altitudes 
1000 km)t This Joule heating is not sufficient to drive the bulk 0*^ to 
escape temperature (*^ 10 eV)* 

2 

2* The thermal energy gain scales as m EJ where m is the mass of the 
dominant F-reglon ion. ITius, the mechanism prefers 0*^ heating. 

3. An initially isotropic Maxwellian distribution is driven to an 
isotropic asymptotic state with a higher temperature by an external Ej^ Im- 
posed on the system. In the transition state, the distribution is non- 
Maxwellian with a superthermal tail. The ion distribution in this period (in 
the vicinity of the t x ^ drift energy y m u^) is anisotropic with respect to 
the magnetic field. 

4. The thermal energy Increase of the entire ion population in the 
transition state is also anisotropic with perpendicular thermal energy aver- 
aging several times that of the parallel thermal energy. This anisotropy 
disappears after several collision times minutes), when the transition from 
initial to asymptotic states is complete. 
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Recently Lockwood axis\ Tithorldge (1981) discussed, from the standpoint of 
observed 0’^/H+ transition altitudes, the necessity to assume an ionospheric O'** 
heating mechanism to raise the 0*** temperature to several eV in the auroral 
region. If Mozer's observation [1981) of large Ej^ in the topside ionosphere 
is confirmed, we believe this can be accomplished by Joule heating in the 
auroral region. We emphasize that our discussion of Joule heating (heating of 
the bulk 0'*' population) is in qualitative agreement with the calculation of 
Redder and Banks (1972) if we apply our results to the large-scale convection 
electric field which maps to tens of mV/ra in the ionosphere, rather than to 
the auroral electric field. Auroral electric fields of larger magnitude ate 
more localized and do not map through to the lower ionosphere [Rich et al., 
1901). A probable reason for this may be the effect of the ionospheric con- 
ductivity profile on the mapping of electric fields inside the ionosphere 
[e.g. Chiu, 19741 • Joule dissipation of strong electric fields in the topside 
ionosphere, such as proposed here, may be a second reason for the absence of 
strong electric fields in the ionosphere. 

Optical observations of the auroral 6300 A line of 0(^D) indicate that 

the intensity is too high (by about one order of magnitude) to be explained by 

either electron Impact on atomic oxygen and/or by dissociative recombination 
“4" .4* 

of O 2 [Sharp et al., 1979]. If the bulk of the 0“ ions on the topside can be 
heated to > 2 eV or if there is sufficient superthermal flux of 0'*’, Increased 
population of the 0(^D) state can easily be accomplished by the charge ex- 
change interaction between the hot (> 2 eV) O'*" and cold atomic oxygen - the 
reaction-product atomic oxygen can easily be in the state [private communi- 
cation, A. B. Christensen] • 

Possibly our model for producing superthermal 0'*’ by auroral electric 
fields is a first step toward a theory of auroral ion beams and conics. While 
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wo must conatontly bo roinlnded of the extcorae slmpUficatioti of the model, the 
Bupotthormal pact of the kinetlc-rooponBO dlotrlbutiotv function of our model 
PQ 8 BO 8 S 0 S ©SBontially the attributes of the conic distribution since the 
particles that expand upwards alons the magnetic field will acquire a 
parallel velocity from the conservation of the magnetic moment. It remains to 
be shown, however, that the model supertherraal flux at keV ouergles Is 
sufficient to account for the observed flux. We hope to do a thorough numeri- 
cal analysis of the distribution in the future. Clearly our model is not 
intended to deal with high altitude (» 1000 Icm) phenomena where, we expect, 
little collisional influence and wave-particle interaction to be important. 
Observationally, it would be interesting to see if measurements of low-energy 
anisotropic superthermal populations of ions at low altitudes (~ 400 km) 
[Whalen et al., 19781 can be extended to energies as low as several aV, as in 
Klurapar [1979]. We believe, as do liockwood and Titheridge [1981], that a 
complete description of auroral 0*^ energization and injection into the magnet- 
osphere must begin with kinetic processes in the F-region ionosphere itself. 

After completion of this study, it was brought to our attention that 
similar ideas and formulation were discussed in a preprint issued by the Space 
Research Institute of the Soviet Academy of Sciences [Zakharov et al.,1980]. 
The discussions of the Russian work do not involve solutions of the Krook 
model os is* done here. 
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APPENDIX A; TRANSFORMATION OF THE BOLTZMANN OPERATOR 

Application of the transformation (x,Vjf;) (13) - (16), to 
the components of the Boltzmann operator V t) * ^^^® brackets 
indicate functionals, yields the sought -for identity 


^{I.v.t) - It f(J'.v'.t) 

provided the force is given by 

f (v) ■ e + e(v X Bz)/c 


(A-1) 


(A-2) 


as assumed in oUr model* 

The calculation is straightforward but somewhat tedious, so we shall only 
provide the salient points here. Direct differentiation gives 


(x,v,t) 

dt 


3f(x’.v'.t) ^ ax'. 



V 


f(x',v',t) 



f(x',^,t)(A-3) 


where V and V^, are gradient operators with respect to x' and v' respective- 
ly. From (13) - (16), we have 


3x' 

3t 


- v(x’,v',t) 


(A-4) 


- t(J') 


(A-5) 


Next , 
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bOCttUBO 


o 


V 


Vf(x,v,t) 


v(x' |V' ,t) 


V’f(x',vSt) 


(A"6) 


V* • 3x /Dx’ - V' . Application of (13) “ (16) once mote gives 
1 111 


t(?) . £(I,v,t) 


t(J') 


V^, £(5',?',0 


(A-7) 


Ilie Identity (A-l) Is obtained by using (A-3) - (A-7). 


APPENDIX B! CONSERVATION OP PARTICLES 


It l8 pointed out in Section III thoU the choice of model initial and 
asymptotic distribution functions, (12) and (7) respectively, must be self- 
consistent with the imposed model constraint ^ which amounts to (21) 
via (4); 

n(x',t) / d^v' f(x’,v*,t) “ Oq (b"!) 

where f is given in terms of the forms (7) and (12) by (20). Here we prove 
(B-1) by direct integration. 

Subatltutlon of (20) Into (B-1) yields two integrations over d^v’ . The 
velocity space Integration over f (x’ ,v’ ,0), (12), yields straightforwardly 
no« The velocity space Integration over Iq looks very complicated but inspec- 
tion of (24) shows that 


(v - u)^ ■ [v* - w(t)]^ (B-2) 

where w(t) does not depend on v'. The Maxwellian form of Iq now allows 
/ d v’ Iq to be Integrated, yielding a result independent of w. Thus, (B-1) 
yields 


n(x%t) « Oq e + Oq V e '*'■ / dt c “ »q 


VT 
e “ n. 


(B-3) 


preceding page blank not filmed 
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